www.gusucode.com > 小波分析理论与Matlab 7实现 > 小波分析理论与Matlab 7实现/code/Ch17/例程17-2.m

    %装载语音信号
N=1024;
s=wavread('m?l.wav',N);
figure(1);
plot(1:N,s,'LineWidth',2);
xlabel('时间 n');
ylabel('幅值 A');

%选用小波函数
wavelet='db5';
%分解级数
level=5;
%压缩参数
alpha=1.5;
%阈值类型
sorh='h';
%小波分解
[c,l]=wavedec(s,level,wavelet);
%计算阈值
[thr,nkeep]=wdcbm(c,l,alpha);
%使用硬阈值压缩信号
[xc,cxc,lxc,perf0,perfl2]=wdencmp('lvd',c,l,wavelet,level,thr,sorh);
%对系数编码
cC=encode(cxc);
%传输的系数
tC=cC;
%传输的系数向量长度
tL=length(tC);
%零系数的百分比
pzeros=perf0
%保留的能量大小
pnormen=perfl2
%编码的压缩比
compratio=length(s)/length(tC)